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Abstract 



Massless quark and antiquark jet production at large transverse momen- 
tum in the coherent diffractive processes at hadron colliders is calculated in 
the two-gluon exchange parametrization of the Pomeron model. We use the 
helicity amplitude method to calculate the cross section formula. We find 
that for the light quark jet production the diffractive process is related to 
the differential off-diagonal gluon distribution function in the proton. We 
estimate the production rate for this process at the Fermilab Tevatron by ap- 
proximating the off-diagonal gluon distribution function by the usual diagonal 
gluon distribution in the proton. And we find that the cross sections for the 
diffractive light quark jet production and the charm quark jet production are 
in the same order of magnitude. We also use the helicity amplitude method 
to calculate the diffractive charm jet production at hadron colliders, by which 
we reproduce the leading logarithmic approximation result of this process we 
previously calculated. 
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I. INTRODUCTION 



In recent years, there has been a renaissance of interest in diffractive scattering. These 
diffractive processes are described by the Regge theory in terms of the Pomeron (IP) exchange 
The Pomeron carries quantum numbers of the vacuum, so it is a colorless entity in QCD 
language, which may lead to the "rapidity gap" events in experiments. However, the nature 
of Pomeron and its reaction with hadrons remain a mystery. For a long time it had been 
understood that the dynamics of the "soft pomeron" is deeply tied to confinement. However, 
it has been realized now that how much can be learned about QCD from the wide variety 
of small- x and hard diffractive processes, which are now under study experimentally. In 
Refs. the diffractive J/ip and T production cross section have been formulated in 

photoproduction processes and in DIS processes in perturbative QCD. In the framework 
of perturbative QCD the Pomeron is represented by a pair of gluon in the color-singlet 
sate. This two-gluon exchange model can successfully describe the experimental results 
from HERA §. 

On the other hand, as we know that there exist nonfactorization effects in the hard 
diffractive processes at hadron colliders First, there is the so-called spectator effect 

||, which can change the probability of the diffractive hadron emerging from collisions intact. 
Practically, a suppression factor (or survive factor) u Sp n is used to describe this effect [ID 



Obviously, this suppression factor can not be calculated in perturbative QCD, which is now 
viewed as a nonperturbative parameter. Typically, the suppression factor Sp is determined 
to be about 0.1 at the energy scale of the Fermilab Tevatron |J. Another nonfactorization 
effect discussed in literature is associated with the coherent diffractive processes at hadron 
colliders 0, in which the whole Pomeron is induced in the hard scattering. It is proved in 
that the existence of the leading twist coherent diffractive processes is associated with a 
breakdown of the QCD factorization theorem. 

Based on the success of the two-gluon exchange parametrization of the Pomeron model 
in the description of the diffractive photoproduction processes at ep colliders J2],|3|,|5j , we may 
extend the applications of this model to calculate the diffractive processes at hadron colliders 
in perturbative QCD. Under this context, the Pomeron represented by a color-singlet two- 
gluon system emits from one hadron and interacts with another hadron in hard process, in 
which the two gluons are both involved (as shown in Fig. 1). The partonic process is plotted 
in Fig.2, where there are nine diagrams in the leading order of perturbative QCD. Therefore, 
these processes calculated in the two-gluon exchange model are just belong to the coherent 
diffractive processes in hadron collisions. Another important feature of the calculations 
of the diffractive processes in this model recently demonstrated is the sensitivity to the 
off-diagonal parton distribution function in the proton [|TT|]. 

Using this two-gluon exchange model, we have calculated the diffractive J /if) production 
12 1, charm jet production []T3], massive muon pair and W boson productions in hadron 



collisions. These calculations show that we can explore much low x (off-diagonal) gluon dis- 
tribution function and study the coherent diffractive processes at hadron colliders through 
these processes. In this paper, we will calculate the light quark (massless) jet production 
at large transverse momentum in the coherent diffractive processes at hadron colliders by 



using the two-gluon exchange model. In the calculations of Refs. [|12]-[1J], there always is a 
large mass scale associated with the production process. That is for J/ip production, 



2 



m c for the charm jet production, M 2 for the massive muon production (M 2 is the invariant 
mass of the muon pair) and Mjy for W boson production. However, in the light quark jet 
production process, there is no large mass scale. So, for the light quark jet production, 
the large transverse momentum is needed to guarantee the application of the perturbative 
QCD. Furthermore, the experience on the calculations of the diffractive di-quark photopro- 



duction JT^] shows that the light quark jet production in the two-gluon exchange model has 
a distinctive feature that there is no contribution from the small Ij, region {Ij, < kj<) in the 
integration of the amplitude over Ij,. So, the expansion (in terms of l\jM\) method used 
in Refs. can not be applied to the calculations of light quark jet production. In the 

following calculations, we will employ the helicity amplitude method to calculate the am- 
plitude of the diffractive light quark jet production in hadron collisions. We will show that 
the production cross section is related to the differential (off-diagonal) gluon distribution 



function in the proton as that in the diffractive di-quark jet photoprodution process |L5 
(On the other hand, we note that the cross sections of the processes calculated in Refs. 
|I2H14| are related to the integrated gluon distribution function in the proton). 

The diffractive production of heavy quark jet at hadron colliders has also been studied by 
using the two-gluon exchange model in Ref. ||16|| . However, their calculation method is very 
different from ours []. In their calculations, they separated their diagrams into two parts, 
and called one part the coherent diffractive contribution to the heavy quark production. 
However, this separation can not guarantee the gauge invariance fll3f| . In our approach, 
we follow the definition of Ref. 0, i.e., we call the process in which the whole Pomeron 
participants in the hard scattering process as the coherent diffractive process. Under this 
definition, all of the diagrams plotted in Fig. 2 for the partonic process gp — > qqp contribute 
to the coherent diffractive production. 

The rest of the paper is organized as follows. In Sec. II, we will give the cross section for- 
mula for the partonic process gp —>■ qqp in the leading order of perturbative QCD, where we 
employ the helicity amplitude method to calculate the amplitude for this process. In Sec. Ill, 
we use this helicity amplitude method to recalculate the diffractive charm jet production 
process gp — > ccp, by which we will reproduce the leading logarithmic approximation result 
for this process previously calculated in Ref. |T3]| . In Sec. IV, we will estimate the production 
rate of diffractive light quark jet at the Fermilab Tevatron by approximating the off-diagonal 
gluon distribution function by the usual diagonal gluon distribution function in the proton. 
And the conclusions will be given in Sec.V. 

II. THE CROSS SECTION FORMULA FOR THE PARTONIC PROCESS 

For the partonic process gp — > qqp, in the leading order of perturbative QCD, there are 
nine diagrams shown in Fig. 2. The two-gluon system coupled to the proton (antiproton) in 
Fig. 2 is in a color-singlet state, which characterizes the diffractive processes in perturbative 
QCD. Due to the positive signature of these diagrams (color-singlet exchange), we know 
that the real part of the amplitude cancels out in the leading logarithmic approximation. 



1 For detailed discussions and comments, please see [13] 
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To get the imaginary part of the amplitude, we must calculate the discontinuity represented 
by the crosses in each diagram of Fig. 2. 

The first four diagrams of Fig. 2 are the same as those calculated in the diffractive photo- 
production processes. But, due to the existence of gluon-gluon interaction vertex in QCD, 
in the partonic process gp —>■ qqp, there are additional five diagrams (Fig.2(5)-(9)). These 
five diagrams are needed for complete calculations in this order of QCD. 

In our calculations, we express the formulas in terms of the Sudakov variables. That is, 
every four-momenta ki are decomposed as, 

hi = a { q + fap + k iT , (1) 

where q and p are the momenta of the incident gluon and the proton, q 2 = 0, p 2 = 0, and 
2p ■ q = W 2 = s. Here s is the cm. energy of the gluon-proton system, i.e., the invariant 
mass of the partonic process gp — > qqp. ai and fa are the momentum fractions of q and p 
respectively, k^ is the transverse momentum, which satisfies 

k iT ■ q = 0, k iT ■ p = 0. (2) 

All of the Sudakov variables for every momentum are determined by using the on-shell 
conditions of the momenta represented by the external lines and the crossed lines in the 
diagram. The calculations of these Sudakov variables are similar to those in the diffractive 
charm jet production process gp — *> cap. And all of these Sudakov variables for the process 
gp — > qqp of the light quark jet production can take their values from the corresponding 
formulas in Ref. |H| after taking m c — > 0. For convenience, we list all of the Sudakov 
variables for the diffractive process gp — > qqp in the following. 

For the momentum u, we have 

Ml o 

at u = 0, fa = xip = , u T = t = 0, (3) 

s 

where M\ is the invariant mass squared of the diffractive final state including the light 
quark and antiquark jets. For the high energy diffractive process, we know that M 2 <C s, so 
we have fa (xjp) as a small parameter. For the momentum k, 

k 2 

a k (l + a k ) = - — fa = -a k /3 u , (4) 
M x 

where k? is the transverse momentum of the out going quark jet. For the loop mentum Z, 

Oil = , 

s 

n 2(k T , l T ) ~ l 2 T r 

fa = , for Diag.l, 3, 5, 

a k s 

n 



2(k T J T ) + l 2 T 
(l + a k )s 
M\-l 2 T 



for Diag.2, 4, 6, 

for Diag.7, 8, 9. (5) 
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Using these Sudakov variables, we can give the cross section formula for the partonic 
process gp — »■ qqp as, 

da(gp^qqp) dM 2 x d 2 k T da k fcf, 2 

It lt =° = Wftr3^ i(a ' (1 + ° fc) + M| ) E ^ ' (6) 

where A is the amplitude of the process gp — > ggp. We know that the real part of the 
amplitude ^4 is zero, and the imaginary part of the amplitude A(gp —>■ qqp) for each diagram 
of Fig. 2 has the following general form, 

ImA = C F (T$ [ x ui(k + q)T,v 3 {u - k), (7) 

where Cp is the color factor for each diagram, and is the same as that for the gp — > cap 
process a is the color index of the incident gluon. T ^ represents some 7 matrices 



including one propagator. F in the integral is defined as 

3 



where 



^-/jix'y-ii), (8) 



J[x,x,l T )- , (Jj 



where the function G(x', x"\ k\) is the so-called off-diagonal gluon distribution function [11 
Here, x' and x" are the momentum fractions of the proton carried by the two gluons. It is 
expected that at small x, there is no big difference between the off-diagonal and the usual 



diagonal gluon densities [17]]. So, in the following calculations, we estimate the production 
rate by approximating the off-diagonal gluon density by the usual diagonal gluon density, 
G(x' ,x";Q 2 ) w xg(x,Q 2 ), where x = x F = M\js. 

In Ref. we calculate the amplitude (0) for the heavy quark diffractive production 
processes by expanding T M in terms of l\. However, in the light quark jet production 
process gp — > qqp, the expansion method is not yet valid. According to the result of [TB|], the 
production cross section is proportional to the heavy quark mass. If we apply this formula 
to the light quark jet production, the cross section will be zero. That is to say, the expansion 
of the amplitude in terms of Z^, in which the large logarithmic contribution comes from the 
region of l\ <C M x , is not further suitable for the calculation of the cross section for massless 
light quark jet production. Furthermore, the experience of the calculation of the diffractive 
light quark jet photoproduction process •yp — > qqp (for Q 2 = 0) |L5| indicates that there is no 
contribution from the region of l\ < k\ in the integration of the amplitude over l\. In the 
hadroproduction process gp — > qqp, the situation is the same. So, the expansion method, 
in which l\ is taken as a small parameter, is not valid for the calculations of the massless 
quark production processes. 



In the following, we employ the helicity amplitude method [0] to calculate the amplitude 
Eq.(^). For the massless quark spinors, we define 

ti±(p) = -^(l±7fe)«(p)- (10) 



5 



For the polarization vector of the incident gluon, which is transversely polarized, we choose, 



e± = -^(0,l,±z,0). (11) 

The helicity amplitudes for the processes in which the polarized Dirac particles are involved 
have the following general forms [18] , 



u±{p f )Qv T { Pl ) = J — (12) 

A \J{ n -Pi)( n m Pf) 

where n is an arbitrary massless 4-vector, which is set to be n = p in the following calcula- 
tions. Using this formula (|12D, the calculations of the helicity amplitude A(X(g), X(q), X(q)) 
for the diffractive process gp —>■ qqp is straightforward. Here A are the corresponding he- 
licities of the external gluon, quark and antiquark. In our calculations, we only take the 
leading order contribution, and neglect the higher order contributions which are proportional 

M 2 

to f3 u = because in the high energy diffractive processes we have p u -C 1. 

For the first four diagrams, to sum up together, the imaginary part of the amplitude 
A(±,+,-) is 

t 41234/, \ a/1 \ \ r f d? l T £ , , „ fl w 2e^> -k T 1 e& ■ (k T - l T ) 
Inu4 1234 (±,+,-) =a 2 k {l + a k )Nx / -—f{ x >, x "-l 2 T )(--— ra ^ tT y r ^ J ), 

(13) 

where | and — ^ are the color factors for Diag.1,4 and Diag.2,3 respectively, and M is defined 

as 

, * y s T? r (14) 

-a k (l + a k ) 

The other helicity amplitudes for the first four diagrams have the similar form as (|i~3|). 
The amplitude expression Eq. (|13|) is the same as that for the photoproduction process 
■yp — > qqp previously calculated in Refs. [|4|,|T5f except the difference on the color factors. 
In the diffractive photoproduction process, the color factors of these four diagrams are the 
same (they are all |), while in hadroproduction process the color factors are no longer 
the same for these four diagrams. It is instructive to see what is the consequence of this 
difference. We know that the amplitude of the diffractive process in Eq. ([/]) must be zero 
in the limit l\ —>■ 0. Otherwise, this will lead to a linear singularity when we perform the 
integration of the amplitude over l\ due to existence of the factor \j{\\) 2 in the integral of 
Eq. (0) Hl3| . This linear singularity is not proper in QCD calculations. So, we must first 



exam the amplitude behavior under the limit of l\ — >• for all the diffractive processes in 
the calculations using the two-gluon exchange model. From Eq. (0), we can see that the 
amplitude for the diffractive photoproduction of di-quark jet process is exact zero at l\ -> 0. 
However, for the diffractive hadroproduction process gp — > qqp the amplitude for the first 
four diagrams is not exact zero in the limit l\ — > due to the inequality of the color factors 
between them. So, for gp — > qqp process there must be other diagrams in this order of 
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perturbative QCD calculation to cancel out the linear singularity which rises from the first 
four diagrams. The last five diagrams of Fig. 2 are just for this purpose. 
For example, the contributions from Diags.5 and 8 are 

IhU-(±, +, -) = ««1 + a k )N x / ||/(*<, *»; + 

I'tJ 4 (kt - (i + 0L k yi T y 

(15) 

and the contributions from Diags.6 and 9 are 

W-(±,+,-) =a2(l + a^x / ^ / (^^;4)(^ ^ ) r fe - a f ) ), (16) 
and the contribution from Diag.7 is 

W 7 (±,+,-)=^(l + a fe )Arx J ^f(jy ; 1*)(-—*.). (17) 



From the above results, we can see that the contributions from Diags.5-9 just cancel out the 
linear singularity which rises from the first four diagrams. Their total sum from the nine 
diagrams of Fig. 2 is free of linear singularity now. 

Finally, by adding up all of the nine diagrams of Fig. 2, the imaginary parts of the 
amplitudes for the following helicity sets are, 

lmA(±, +, -) = a 2 k (l + a k )M x T (±) , 

Invl(±, -, +) = a k (l + a k ) 2 Af x T (±) , (18) 



where 



-J (P T y nx > x '^h 9 j k\ 36 (k T -i T y 

1 + a k ef^ ■ (k T - (1 + OL k )l T ) a k ■ (k T - a k h) , Mn s 

+ ~ 7? rT^J- l 1J ) 



4 {k T - (i + ak )l T y 4 (k T -a k l T) 

And the amplitudes for the other helicity sets are zero in the light quark jet production 
process. From the above results, we can see that in the integration of the amplitude the linear 
singularity from different diagrams are canceled out by each other, which will guarantee there 
is no linear singularity in the total sum. 

Another feature of the above results for the amplitudes is the relation to the differential 
off-diagonal gluon distribution function f(x', x"; To get the cross section for the diffrac- 
tive process, we must perform the integration of Eq. (0). However, as mentioned above that 
there is no big difference between the off-diagonal gluon distribution function and the usual 
gluon distribution at small x, so we can simplify the integration of (|T^p by approximating 
the differential off-diagonal gluon distribution function f (x' , x" ; l^) by the usual diagonal 
differential gluon distribution function f g (x; 

After integrating over the azimuth angle of It, the integration will then be 
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TW-xf^Z, (20) 



where 



X ~~ / 7r^J9\ x i l T)Y^Ko 7TTT2 727 > + 



(i2)2^v-'-iVL 36 v 2 2 |fcf,-Z£| / ' 4 v 2 2|^-(l + a fc )4| 



•)]■ (21) 



4^2 2\k\-a k l\ 

Comparing the above results with those of the photoproduction process •yp — > qqp fl4yi5f. 
we find that the amplitude formula for the diffractive light quark jet hadroproduction pro- 
cess gp — > qqp is much more complicated. However, the basic structure of the amplitude, 
especially the expression for the integration X is similar to that for the photoproduction 
process. In the integration of fl2"T|), if I 2 < k\ the first term of the integration over l\ will 
be zero; if I 2 < fcf-/(l + ak) 2 the second term will be zero; if I 2 < k\ja\ the third term will 
be zero. So, the dominant regions contributing to the three integration terms are I 2 ~ k\, 
l\ ~ k\j{\ + ak) 2 , and l\ ~ k\ja\ respectively. Approximately, by ignoring some evolu- 
tion effects of the differential gluon distribution function f g {x\l\) in the above dominant 
integration regions, we get the following results for the integration X, 

I = i 4) + {1 ~^f:M: **) - 4u.r: §)]. (22) 



A' 



Obtained the formula for the integration X, the amplitude squared for the partonic 
process gp — > qqp will be reduced to, after averaging over the spin and color degrees of 
freedom, 



]A\ 2 = ^(4vr) WE^(1 - M). (23) 

And the cross section for the partonic process gp — > qqp is 

rfa( OT -> qqp) _ r a ^ 9^tt 2 1 2A ^| T ,2 , ?4 n 



The integral bound M\ > 4k 2 -, above shows that the dominant contribution of the integration 
over M\ comes from the region of M\ ~ 4k 2 -,. Using Eq. @, this indicates that in this 
dominant region ak is of order of 1. So, in the integration X the differential gluon distribution 
function f g (x; Q 2 ) of the three terms can approximately take their values at the same scale 
of Q 2 = k\. That is, the integration X is then simplified to 

J ~ 36Mi ^**)- (25) 

Numerical calculations show that there is little difference (within 10% for kx > 5 GeV) 
between the cross sections by using these two different parametrizations of X, Eq. ( |2"2"| ) and 
Eq. (p5|) . So, in Sec. IV, we use Eqs. ( p4|) and ( p5|) to estimate the diffractive light quark jet 
production rate at the Fermilab Tevatron. 
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III. RECALCULATE THE HEAVY QUARK JET PRODUCTION USING THE 

HELICITY AMPLITUDE METHOD 



For a crossing check, in this section we will recalculate the diffractive heavy quark jet 
production at hadron colliders by using the helicity amplitude method. In Ref. |0|, we 
have calculated this process in the leading logarithmic approximation of QCD, where we 
expanded the amplitude in terms of l\. Now, if we use the helicity amplitude method, we 
donot need to use the expansion method for the Y ^ factor in Eq. (^) as in ||13fl . We can 
firstly calculate the amplitude explicitly by using the helicity amplitude method. 

However, for the massive fermion, the amplitude formula is more complicated. Following 
Ref. [fL9| , we first define the basic spinors u±(k ) as, 

u + (k ) =}6iU-(ko), 
u_{k )u_{k ) = ^(l-Ts) #0, (26) 

where the momenta k and k% satisfy the followoing relations, 

k -k = 0, kx-ki = -1, k -h = 0. (27) 
Using Eqs. Q2"S| ) and fl2"T|), we can easily find that the spinor u + (k ) satisfies 

u + (k )u+(k ) = -(l+ l5 ) ft. (28) 



Provided the basic spinors, we then express any spinors u(pi) in terms of the basic ones, [19 

= W + >"<>M^ (29) 

It is easily checked that these spinors satisfy Dirac's equations. Now, for the massive 
fermions, the helicity amplitudes for the processes involving Dirac particles have the fol- 
lowing general forms, 



u+{Pf)Qv-{pi 
U-(Pf)Qv+(pi 
u + (p f )Qv + (pi 



Tr[Q(tfj - mj)(l - 75) fojjtf + m f ) 

Tr\Q(jk - + 7s) jjo(|j/ + m f ) 

Tr\Q(Ji - mj)(l - 75) fa fajjfo + m f ) 
^(n ■ Pi)(n • p f ) 



n-{Pf)Qv-{pd = TrMA-^Ml-riWf + rnf)^ (3Q) 

4^ (n ■ Pi)(n ■ p f ) 



where m 8 and my are the masses for the momenta Pi and pf respectively, where pf = mf , p 2 j 
m 2 f. From the above equations, we can see that for the massless fermions (m^ = rrif = 1 
the formula Eq. fl30|) will then turn back to the formula Eq. (|T2|). 
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To calculate the imaginary part of the amplitude Eq. (0) for the partonic process gp — > 
ccp, a convenient choice for the momenta ko and ki is, 

k =p, ki = e, (31) 

where the vector e is the polarization vector for the incident gluon defined in Eq. (|Tl~l). 
Using the formula Eq. ([!(]) and the above choice for the momenta ko and k\, the helicity 
amplitudes for Eq. (|7|) will then be, 

lmA(±, +, -) = a 2 k (l + a k )N x 

W(±, -, +) = a k {l + a k ) 2 M x T C W, (32) 

7T777 

lmA(±, +, +) = W(±, -, -) = a k (l + a k )M x (33) 
where J\f is the same as in Eq. (TTJ]), and the integrations 7^ an d are defined as 

(±) _ r^r , „ 2 1 2 j^-jr 1 ■ (fcr ~ h) 
c J (l 2 T r [X,X ,lT)[[ 2 + 36 m 2 + ( ^_f T )2 

1 + a k e (± ) • (A; T - (1 + o^t) + a fc e^j ■ (fc T - o^r) ^ ^ 



4 m 2 + (& T _ (i + « fc )/ T )2 4 m 2 + (fcr _ afc/r)2 J 

J2L 1 O 1 J ]_ 



X ' c ~ nl (l^ f{x '' X "' j2T)[{ 2 9 ] k 2 + 



+ ^ ,J - (35) 



4 m 2 + (& T _ (i + ak )i T )2 4 m 2 + (A;t _ afc /t)2 j 

If we approximate the differential off-diagonal gluon distribution function f(x', x"\ l 2 ^) by the 
usual diagonal differential gluon distribution function f g (x;l T ), the above integrations will 
then be reduced to, after integrating over the azimuth angle of It, 



where 



T c (±) = 7re< ±) • k T X c , (36) 



d 2 l T 2 5 1 5 1 1 k 2 T -m 2 c -l 2 T 
JjffihW l T)lj^^T ~ 35 ^ ~ 72 ~\ 

1 + a k fc| - mg - (1 + a fc ) 2 ^ a fc k\ - m 2 c - a\l 2 T ^ 



m 2 8 m% 



where 



m 2 T = k 2 T + m 2 c , m\ = yj \m\ + l\) 2 - k\l\, m\ = \] {m 2 t + (1 + a k ) 2 l T ) 2 - (1 + a k ) 2 k 2 T l 2 T , 
m\ = ^(m 2 + a 2 l 2 T ) 2 -a 2 k 2 l 2 , (38) 

and 

,_ fd% 2 5 1 1 1 l + a k 1 OfcJ_ 

C "i (Z2)2^^'^L 18m 2 36m 2 4 m 2+ 4m 2J- ^ 
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So, the amplitude squared for the partonic process gp — > ccp will then be reduced to, after 
averaging over the spin and color degrees of freedom, 

R 2 = lo? a (An) W^[(l - |J )k 2 T \l c \ 2 + ml\l'f]. (40) 

From the above results, we can see that the integrals of Eqs. fl37D and ( |3~9f ) are proportional 
to in the limit of l\ —>■ 0. That is to say that there exist large logarithmic contributions 



from the integration region of 1/ i% <C l\ <C for the integration over l\ as in Ref . |T3 
So, we can expand the integrals of Eqs. fl3"T|) and (|3D| ) in terms of l\ to get the leading 
logarithmic contribution to the amplitude. In the limit of 1^ — > 0, the parameters m 2 , m 2 
and m\ scale as, 

1 ^\_ [x _ml-k\ l 2 T] 



m 2 m\ m\ m\ 



1 m 2 - 4 (1 + cxfe) 2 ^ 



2 2 L 2 2 J ' 

— ^ — 1- c 2 r fc 2 T . 41 



Under this approximation, the integrations Eqs. (£37]) and ([39]) will then be related to the 
integrated gluon distribution function xg[x\ Q 2 ), 

2m 2 

J ' c W 3 6 r4") 3 Ml (10Mi " 27m ^ (x; m ^ } - (43) 



Substituting the above results into Eq. (f40|) , we can then reproduce the leading logarithmic 
approximation result for the diffractive charm jet production process at hadron colliders 



which has been calculated in [13]. 



IV. NUMERICAL RESULTS 

Provided with the cross section formula (|23|) for the partonic process gp —>■ qqp, we can 
calculate the cross section of the diffractive light quark jet production at the hadron level. 
However, as mentioned above, there exists nonfactorization effect caused by the spectator 
interactions in the hard diffractive processes in hadron collisions. Here, we use a suppression 
factor Ts to describe this nonfactorization effect in the hard diffractive processes at hadron 
colliders [§]. At the Tevatron, the value of Ts may be as small as Ts ~ 0.1 P||. That is to 
say, the total cross section of the diffractive processes at the Tevatron may be reduced down 
by an order of magnitude due to this nonfactorization effect. In the following numerical 
calculations, we adopt this suppression factor value to evaluate the diffractive production 
rate of light quark jet at the Fermilab Tevatron. 

The numerical results of the diffractive light quark jet production at the Fermilab Teva- 
tron are plotted in Fig.3 and Fig.4. In our calculations, the scales for the parton distribution 
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functions and the running coupling constant are both set to be Q 2 = k\. For the parton 
distribution functions, we choose the GRV NLO set [Q. In Fig.3, we plot the differential 
cross section da/dt\ t=Q as a function of the lower bound of the transverse momentum of the 
light quark jet, /crmm- This figure shows that the cross section is sensitive to the transverse 
momentum cut A/rmm- The differential cross section decreases over four orders of magnitude 
as /cTmin increases from 5 GeV to 15 GeV. 

It is interesting to compare the cross section of the diffractive light quark jet production 
with that of the diffractive charm quark jet production |fl3| , which is also shown in Fig.3. The 



two curves in this figure show that the production rates of the light quark jet and the charm 
quark jet in the diffractive processes are in the same order of magnitude. However, we know 
that the cross section of diffractive heavy quark jet production is related to the integrated 
gluon distribution function, while the cross section of the light quark jet production is related 
to the differential gluon distribution function. 

In Fig. 4, we plot the differential cross section da/dt\ t=0 as a function of the lower bound 
of the momentum fraction of the proton carried by the incident gluon Xi m ; n , where we set 
^Tmin = 5 GeV. This figure shows that the dominant contribution comes from the region of 
x\ ~ 10~ 2 — 10 _1 . This property is the same as that of the diffractive charm jet production 
at the Tevatron UTS . 



V. CONCLUSIONS 

In this paper, we have calculated the diffractive light quark jet production at hadron 
colliders in perturbative QCD by using the two-gluon exchange model. We find that the 
production cross section is related to the squared of the differential gluon distribution func- 
tion dG(x;Q 2 )/dlnQ 2 at the scale of Q 2 ~ k^, where kx is the transverse momentum of 
the final state quark jet. For a crossing check, we also used the helicity amplitude method 
to calculate the diffractive charm jet production process at hadron colliders, by which we 
could reproduce the leading logarithmic approximation result for this process previously 



calculated in Ref. [13). We have also compared the production rate of the light quark jet 
in the diffractive processes with that of the heavy quark jet production, and found that the 
production rates of these two processes are in the same order of magnitude. 

As we know, the large transverse momentum dijet production in the diffractive processes 
at hadron colliders is important to study the diffractive mechanism and the nature of the 
Pomeron. The CDF collaboration at the Fermilab Tevatron have reported some results 
on this process |J. In this paper, under the two-gluon exchange model in perturbative 
QCD we have calculated the light quark jet production in the diffractive processes. In a 
forthcoming paper, we will calculate the diffractive gluon jet production in hadron collisions 
by using the two-gluon exchange model, which will complete the calculations of the diffractive 
dijet production in perturbative QCD. These results then can be used to compare with the 
experimental measurements to test the valid of the perturbative QCD description of the 
diffractive processes at hadron colliders and may be used to study the diffractive mechanism 
and the factorization broken effects. 
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Figure Captions 



Fig.l. Sketch diagram for the diffractive quark jet production at hadron colliders in pertur- 
bative QCD. 

Fig. 2. The lowest order perturbative QCD diagrams for partonic process gp — > qqp. 

Fig. 3. The differential cross section da/dt\ t =o for the light quark jet production in the 

diffractive processes as a function of fcrmin at the Fermilab Tevatron, where kxmin is the 

lower bound of the transverse momentum of the out going quark jet. For the charm quark 

jet production cross section formula, we take from Ref. and set m c =1.5 GeV. 

Fig. 4. The differential cross section da/dt\ t=0 for the light quark jet production as a function 

of xi m i n , where £i m i n is the lower bound of X\ in the integration of the cross section. 
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